function [mu_T, L_T] = rnd_time_hyperparams(V)
[D,K] = size(V);

beta0 = 1;
nu0 = D;
rh0 = zeros(D,1);
W0 = eye(D);

mu0_ = (V(:,1) + beta0*rh0)/(beta0+1);
beta0_ = beta0 + 1;
nu0_ = nu0 + K; 

dV = diff(V,1,2);

W0_ = inv(inv(W0) + (K-1)*fcov(dV,2) + (beta0/(beta0+1))*(V(:,1)-rh0)*(V(:,1)-rh0)');

W0_ = (W0_+W0_')*0.5;
L_T = wishrnd(W0_, nu0_);
Li_T = inv(beta0_*L_T);
mu_T = rnd_mvn(mu0_, Li_T);
